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Preface 


The  question  of  scattering  an  electromagnetic  signal  from  a 
relativistic  electron  beam  which  is  propagating  through  a gaseous 
medium  is  addressed.  The  dependence  of  the  scattering  cross  section 
on  the  angle  of  incidence  is  studied  and  presented.  This  work  has 
been  supported  by  the  Naval  Se3  Systems  Command. 


1 


NSWC/WOL/TR  76-101 


TABLE  O:^  CONTENTS 

I.  INTRODUCTION  4 

II.  THE  PLASMA  CHARACTERIZATION  6 

III.  THE  INCIDENT  AND  SCATTERED  FIELDS  8 

IV.  BOUNDARY  CONDITIONS  10 

V.  THE  DIFFERENTIAL  AND  TOTAL  SCATTERING  CROSS  SECTIONS  16 

VI.  THE  SHORT  WAVELENGTH  LIMIT  WITH  THE  PLASMA  COLUMN 

STATIONARY 18 

VII.  CONCLUSIONS 2 0 

APPENDICES 

A.  THE  PLANE  WAVE  EXPANSION 2 2 

B.  THE  PLANE  WAVE  TRANSFORMATIONS 2 3 

C.  THE  WAVE  NUMBER  IN  THE  PLASMA 2 7 

D.  THE  TRANSVERSE  FIELD  COMPONENTS  AS  FUNCTIONS  OF  THE 

LONGITUDINAL  FIELDS  29 

E.  THE  SCATTERED  FIELD  COMPONENTS  IN  THE  OBSERVERS  FRAME  32 

REFERENCES 38 


T 


NSWC/WOL/TR  76-101 


ILLUSTRATIONS 

FIGURE 

1.  Plane  Wave  Obliquely  Incident  on  Plasma  Column  34 

2.  The  Measurement  of  Phase  by  Counting  Wave  Crests  ....  35 

3.  The  Invariance  of  Phase 36 

TABLE 

1.  The  Dominant  Expansion  Coefficients  for  the  Scattering 

Cross  Section 37 


3 


NSWC/WOL/TR  76-101 


I.  INTRODUCTION 

The  scattering  of  electromagnetic  waves  from  plasmas  has  been 
studied  by  numerous  researchers  in  recent  years.  Wait-*-  considered 
the  case  of  a plane  wave  obliquely  incident  on  a dielectric  cylinder 
of  infinite  length.  The  incident  magnetic  field  is  transverse  to 
the  cylinder  axis.  Results  for  the  scattered  fields  are  cited  and 
limiting  cases  such  as  the  far  field  approximation,  the  highly  con- 
ducting cylinder  and  the  highly  permeable  cylinder  are  discussed. 
Scattering  of  electromagnetic  waves  from  a plasma  column  at  rest  may 
be  studied  through  an  extension  of  Wait's  work  if  one  simply  charac- 
terizes the  plasma  medium  by  a frequency  dependent  dielectric  con- 
stant e(w).  Many  researchers  have  approached  plasma  scattering 
problems  in  this  fashion.  Shapiro^  has  studied  the  problem  of 
scattering  from  a homogeneous  plasma.  Others 3, 4 have  considered 
inhomogeneous  and  anisotropic  plasma  columns.  Yeh  and  Rusch^  have 
looked  at  the  plasma  sheath  that  is  both  anisotropic  and  radially 
inhomogeneous.  Far  field  patterns  as  well  as  the  back  scattering 
cross  sections  were  numerically  obtained.  Messian  and  Vandenplas® 
added  complexity  to  the  plasma  scattering  problem  by  allowing  the 
column  to  have  an  axial  drift  velocity,  v<<c.  Such  motion  couples 
an  incident  plane  wave,  whose  magnetic  vector  is  parallel  to  the 
column  axis  (an  H wave)  to  a scattered  wave  which  has  a component 
parallel  to  the  cylinder  axis  (an  E wave) . This  phenomena 
naturally  follows  from  the  Lorentz  covariant  form  of  the  Maxwell 
equations  and  is  the  plasma  high  frequency  analogue  of  the 


1.  J.  R.  Wait,  "Scattering  of  a Plane  Wave  from  a Circular  Dielectric 
Cylinder  at  Oblique  Incidence,"  Can.  J.  Phys . 3_3»  189,  1955. 

2.  H.  Shapiro,  Antenna  Lab.  Tech.  Kept.  No.  11,  California  Institute 
of  Technology,  1957. 

3.  F.  A.  Albini  and  R.  G.  Jahn,  "Reflection  and  Transmission  of 
Electromagnetic  Waves  at  Electron  Density  Gradients,"  J.  Appl. 

Phys.  22_,  75,  1961. 

4.  P.  M.  Platzman  and  H.  T.  Ozaki,  "Scattering  of  Electromagnetic 
Waves  from  an  Infinitely  Long  Magnetized  Cylindrical  Plasma,” 

J.  Appl.  Phys.  31,  1597,  1960. 

5.  W.  V.  T.  Rusch  and  C.  Yeh,  "Scattering  by  an  Infinite  Cylinder 
Coated  with  an  Inhomogeneous  and  Anisotropic  Plasma  Sheath,  IEEE 
Transactions  on  Antennas  and  Propagation,  Vol.  AP-15,  _3,  4 52,  1967. 

6.  A.  M.  Messian  and  P.  E.  Vandenplas , "High  Frequency  Effect  Due  to 
the  Axial  Drift  Velocity  of  a Plasma  Column,"  Phys.  Rev.  149, 

131,  1966. 
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7-9 

Roentgsn-Eichenwald  current.  This  is  the  current  generated  by 

pulling  a dielectric  slab  between  the  parallel  plates  of  a capacitor, 
which  establishes  a magneto-static  field  perpendicular  to  both  the 
velocity  of  the  slab  and  the  electric  field  of  the  capacitor.  YehlO 
calculated  the  scattered  fields  emanating  from  a plasma  column  mov- 
ing with  an  arbitrary  velocity  where  the  incident  plane  wave  was 
obliquely  incident.  As  previous  investigators  found,  cross-polarized 
scattered  fields  are  induced  and  for  a particular  angle  of  incidence, 
namely  8 = sin“l(v  /c)  (where  v is  the  velocity  of  the  plasma 

medium,  c is  the  speed  of  lightzin  vacuum).  In  particular,  both  E 
and  H scattered  waves  will  be  excited  when  E or  H waves  are  incident 
on  a plasma  cylinder,  either  obliquely  or  normally. 

Hodjat''""''  considered  the  problem  of  a plane  wave  incident  at  an 
angle  upon  an  electron  beam  confined  by  a glass  envelope,  i.e.  a 
dielectric  material.  He  also  considered  the  separate  cases  of  a wave 
incident  on  a dielectric  cylinder,  a hollow  dielectric  cylinder  and  a 

moving  electron  beam,  respectively.  The  coefficients  of  the  scattered 
waves  were  found.  An  experiment  was  also  conducted  to  measure  the  z 
component  scattered  electric  field  with  a dipole  antenna,  and  the 
cross-polarized  scattered  magnetic  field  with  a loop  antenna.  The 
experiments  were  restricted  to  the  cases  of  the  dielectric  cylinder, 
the  hollow  cylinder,  and  the  electron  beam  surrounded  by  the  die- 
lectric cylinder. 

12 

Kriegsmann  has  studied  the  scattering  of  plane  waves  from  a 
cylindrically  confined  cold  plasma  in  the  geometric  optics  limit,  for 
the  E and  B polarization  cases.  Solutions  to  the  Helmbltz  equation 
for  an  inhomogeneous  medium  are  found  using  the  geometric  theory  of 

7.  W.  C.  Roentgen,  "Uber  die  durch  Bewegung  eines  in  homogenen 
elektrische  Felde  befindlichen  Dielectricums  hervorgerufene 
elektrodynamische  Kraft",  Ann.  Phys.  35,  264,  1888  and 
"Beschreibung  des  Apparates  mit  welchem  die  Versuche  iiber  die 
elektrodynamische  Wirkung  bewegter  Dielectricia  ausgefuhrt 
wurden,"  Ann.  Phys.  40_,  93,  1890. 

8.  A.  Eic’nenwald,  "Uber  die  Magnestische  Wirkungen  bewegter  Korper 
in  elektrostatischen  Felde,"  Ann  Phys.  11 , 421,  1903. 

9.  A.  Sommerfeld,  "Electrodynamics,"  (Academic  Press  Inc.,  New  York, 
New  York,  1952),  pp.  283-235. 

10.  C.  Yeh,  "Scattering  Obliquely  Incident  Microwaves  by  a Moving 
Plasma  Column,"  J.  Appl.  Phys.  40_,  5066,  1969  . 

11.  F.  Hodjat,  "Scattering  of  Electromagnetic  Waves  from  a Moving 
Cylindrical  Electron  3eam, " UCEA  doctoral  thesis,  University 
Microfilms  No.  72-2829,  1971. 

12.  G.  Kriegsmann,  "An  Application  of  the  Method  of  Geometric  Optics 
to  the  Scattering  of  Plane  Electromagnetic  Waves  Off  Cylindri- 
cally Confined  Cold  Plasmas,"  UCLA  doctoral  thesis,  University 
Microfilms  No.  74-29-265,  1974. 
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diffraction  developed  by  Keller  and  his  co-workers.  The  results  of 
this  method  compare  well  with  those  of  exactly  soluble  problems, 
giving  the  same  results  asymptotically.  Such  notions  as  rays, 
amplitudes  and  caustics  are  all  fundamental  features  of  such  an 
approach . 

In  this  paper  we  follow  the  general  approach  of  Yeh,1^  however 
the  physical  parameter  of  interest  is  the  differential  scattering 
cross  section  per  unit  length. 

II . THE  PLASMA  CHARACTERIZATION 

We  consider  a plasma  that  consists  of  electrons  only.  These 
electrons  are  free  to  move  in  response  to  an  applied  field  of  fre- 
quency a) . If  the  plasma  is  rareified  enough*.  then  the  collision 
frequence  „c  between  electrons  may  be  negligible  compared  to  the 
frequency  of  the  applied  field.  According  to  the  Drude  theory-^  the 
conductivity  may  be  given  by 

2 

(II.  1)  a = inQe  /mto 

where  n0  is  the  number  of  electrons  per  unit  volume,  e is  the  elec- 
tron charge,  and  m is  its  mass.  Thus  the  plasma  we  are  considering 
may  be  regarded  as  collisionless. 

Since  no  longitudinal  fields  can  exist  in  a conducting  medium 
in  the  absence  of  an  applied  current  density, 13  we  shall  consider  the 

transverse  wave  field  that  varies  as  exp  (iit-r-iwt)  where  k is  the 
propagation  wave  vector  and  may  be  found  from  the  Maxwell  equations 
to  satisfy  the  equation 

(11. 2)  k2  - (|)2(1  + |£, 

where  n is  the  index  of  refraction  of  the  plasma,  and  £ is  its 
dielectric  constant. 

If  some  non-equilibrium  charge  configuration  exists  in  the  plasma, 
some  of  the  electrons  will  rearrange  themselves  in  an  attempt  to 
screen  out  the  resultant  electric  field.  For  example,  there  may 
exist  an  excess  of  electrons  in  some  planar  y-z  region,  a shortage  in 
another  planar  region  and  a uniform  electric  field  Ex  between  the 
charged  planes.  Any  electron  in  the  interplanar  region  will  experi- 
ence a force 

W 

( 11 . 3)  mx  = -eE 

x 


10.  C.  Yeh,  Op.  Cit. 

13.  J.  D.  Jackson,  "Classical  Electrodynamics,"  (John  Wiley  and 
Sons,  Inc.,  New  York,  1963),  pp.  223-226. 
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where  the  electric  field  is  given  by 
(II. 4)  Ex  = [/e0  = n0ex/e0 

where  £ is  the  surface  charge  density  on  the  positively  charged 
plane,  and  is  the  permittivity  of  free  space.  Substitution  of 
eq  (II. 4)  into  ea  (II. 3)  yields 

..  noe2 

(II.  5)  x + (-2 — ) x = 0 

me0 

This  is  just  a harmonic  oscillator  equation  and  we  may  define  the 
plasma  frequency  Wp  as 

2 2 

(II.  6)  io0  = nQe  /me0 

Physically,  the  electron  experiencing  a force  will  move  until  it 
reaches  the  positively  charged  plane,  however  its  kinetic  energy  will 
carry  it  beyond,  continuing  until  its  energy  is  ultimately  recon- 
verted to  potential  form.  It  will  then  repeat  the  motion  in  a 
periodic  fashion.  This  is  known  as  a plasma  oscillation  where  the 
restoring  force  is  the  Coulomb  interaction  and  the  electron's  mass 
is  the  inertia. 

If  we  substitute  eq  (II. 1)  into  eq  (II. 2)  making  use  of  eq  (II. 6) 
we  obtain 

(11. 8)  k2  ~ (i  _ w 2/u>2 ) 

c2  P 

Since  k = nu/c,  we  find  an  expression  for  the  index  of  refraction 
of  a plasma  given  by 

(11. 9)  n2  ~ 1 - a)  2/w2 

P 

For  the  case  oj>wp  (sometimes  referred  to  as  the  underdense  case)  the 
index  of  refraction  is  real  and  the  electromagnetic  waves  simply 
propagate  through  the  medium.  The  other  situation,  namely  w<wp 
(the  overdense  case)  gives  rise  to  an  imaginary  index  and  electro- 
magnetic waves  will  scatter  and  be  reflected  from  this  plasma 
medium.  In  the  general  situation,  where  the  charged  particle  density 
is  a function  of  position,  the  waves  propagate  through  the  under- 
dense region  until  reaching  some  point  where  the  particle  density 
has  increased  sufficiently  to  become  overdense.  At  this  critical 
surface  the  waves  are  reflected.  There  will  be  some  penetration  of 
this  surface  by  the  fields  but  they  will  quite  rapidly  decay  ex- 
ponentially in  the  interior  of  the  critical  surface. 

It  should  be  noted  that  the  expression  for  the  plasma  frequency 
given,  in  eq  (II. 6)  is  valid  for  a cold  plasma,  i.e.  when  the  elec- 
trons thermal  kinetic  energy  has  a small  effect  on  the  plasma  densi- 
ty. When  this  is  not  the  case  one  must  consider  fluctuations  in 
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the  particle  density  and  a consequent  modification  of  the  expression 
for  the  plasma  frequency.  For  our  purposes  the  cold  plasma  approxi- 
mation will  suffice.  Jackson-1--^  has  considered  the  problem  of  the 
warm  plasma  oscillations  and  derives  appropriate  expressions  for 
the  frequencies  of  these  oscillations  in  terms  of  the  already  defined 
cold  plasma  frequency  and  the  density  fluctuations.  These  expres- 
sions are  appropriately  called  "dispersion  relations"  since  they 
relate  the  oscillation  frequency  m to  its  wave  number  k. 

III.  the  incident  and  scattered  fields 

VJe  shall  consider  a plane  wave  propagating  in  such  a direction 
as  to  be  obliquely  incident  on  a plasma  column  that  is  infinite  in 
extent,  aligned  parallel  to  the  z axis.  Fig.  1 clearly  illustrates 
the  geometry  of  the  initial  conditions.  One  may  characterize  the 
plasma  by  specifying  either  its  conductivity  a (u)  or  its  dielectric 
constant  e (w) . Because  the  plasma  conductivity  will  be  finite,  in 
general,  it  is  more  conventient  to  use  the  dielectric  constant  as 
the  characterization  parameter.  One  may  then  treat  the  problem  of 
the  scattering  by  a plasma  column  in  the  same  fashion  as  the  scat- 
tering from  a dielectric  cylinder. 

We  write  the  expression  for  the  incident  plane  wave  with  respect 
to  the  laboratory  S frame 


(III . i)  E<l}  * Eocoseoexp(-ikoycos9o  + ik0zsin60)  e-1"^ 

H(l)  = 0 
z 

where  E and  w are  the  amplitude  and  frequency  of  the  incoming, 
incident  wave  and 

ko  = u/c  = ^posor£ 

where  we  have  chosen  to  work  in  the  MKSA  system  of  units.  In  a 
reference^ rame  S'  that  is  embedded  in  the  column,  moving  with 
velocity  v = vz  \ z relative  to  the  S frame  the  incident  wave  is 
written : 


(III. 2)  E^  = E'cos9'exp(-iko'cos8 'y'  + ik  sinB'z ')  e 


-iw  't 


= f; 


(-1) nJn (k^r 'cos6 ' ) exp ( in0  ' ) 


n=  -c 


13.  J.  D.  Jackson, 

Op.  Cite,  pp.  223-226. 
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(III. 3)  H x = 0 

2 

The  right  hand  side  of  eq  (III.  2)  is  derived  in  Appendix  A.  The 
wave  parameters  measured  in  the  S'  frame  are  related  to  those 
measured  in  the  S frame.  The  relationships  are 


III. 4) 

Uj  — 

Y = 

8 = 

Yoj  (l-3sinS  ) 

(1-S2)~^ 

V /c 
z 

(III. 5) 

k'cos9'=  k cos6 
o o o 

(III. 6) 

sin6 

= (sin9o-3) / (l-8sin6c) 

(III. 7) 

ii 

' o 

“'(lioeo)%  = Yko(l-Ssineo) 

(III. 8) 

ii 

' o 
w 

y EQ (l-2sin0Q) 

(III. 9) 

fe  " 

E'cosffexp  ( ik^sine  'z  ' - iui 

= Eqcos 6Qexp  ( ikosin9Qz  - iiot) 

Equations  (III. 4)  through  (III. 9)  are  derived  in  appendix  B. 

The  z -component  scattered  fields  may  be  written  as 

°o  . 

(111. 10)  £2'!S)(r',e',o')  = FE'V  (-DV^I^r'cosej'111 

n=-« 

(111. 11)  H'ls)  (r'fe',<m  = FE5>1)ny-(ir)  Hn(1)  ^kor'cose^eJ 


where  the  expansion  coefficients  A^  and  will  be  determined  from 

boundary  conditions.  The  Hankel  function  of  the  first  kind  is  used 
to  represent  outgoing  radiation  and  is  given  by 

(III. 12)  (x)  = Jn(x)  + iNn(x) 

where  J (x)  and  N (x)  are  the  ordinarv  Bessel  and  Neumann  functions 
n n 


9 


respectively . 


We  may  rt 
plasma  column  bv 


2T.t  the  z-comoonent  fields  that  penetrate  into  the 


(111-13) 


30 

E'<p)  (r',e'fc')  = fe^  (- 


1)Vj  (X'r')ein,p 
n n 


(III. 14) 


H'ip)  ( r ' , 9 ' , o ' ) 


n=-j 


n V UQ  / n 


and 


1£1 


2 ~ . k 


(III. 15)  = k ( ~ sin‘0") 


o o 


with  C ' and  D ' determined  from  the  boundary  conditions  and  eq 

n n 

(III. 15)  derived  in  appendix  C. 


IV.  BO  UN DA RY  CONDITIONS 


The  total  tangential  maqnetic  and  electric  field  intensities  are 
continuous  across  the  plasma  - vacuum  interface,  i.e.  the  cviir.der  surface. 

n x (E(l)  + E(S)  ) = n x E(p) 


(IV. 1) 


n x (H(l)  + H(s)  ) = n x H(p) 


where  n is  the  unit  normal  vector  associated  with  the  plasma  surface. 

If  we  use  eqs  (III. 2),  (III. 10)  and  (III. 13)  in  the  first  of 
eq  (IV. 1)  we  get 


(IV. 2) 


> , . f i ) - 

J (k  acose  ) + A Hv  ' (k  acosd  ) - C J ( X a)  = 0 


n o 


n n 


Making  us  of  eqs  (III. 3),  (III. 11)  and  (III. 14)  in  the  second  of 
eq  (IV. 1)  yields 

(IV. 3)  (1)  (koacos6 ')  - Dn«Jn(X'a)  = 0 
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So  far  we  have  two  equations  for  the  four  undetermined  coefficients 

A',  Q' , C'  and  D'.  We  shall  find  two  more  equations,  from  the 
n n'  n n ^ 

Maxwell  equations  which  will  give  us  a complete  set  of  equations  from 

which  these  constants  may  be  determined.  The  Maxwell  equations  are: 


(IV. 4) 

V x E = - 

3B  A 

at  ' 

(IV. 5) 

$ . D = 0 

, 7 x H 

supressing  the  superscripts  "s"  and  ”p"  and  prime  "1"  for 
convenience.  We  shall  utilize  the  constitutive  relations  as  well, 
namely 

(IV. 6)  D = eQE  + P = KeQE 

(IV. 7)  B = uq(H  + M)  = KmUQH 

The  quantities  K and  Km  are  the  relative  permittivity  and  permeability 
respectively.  For  example,  in  vacuum  K and  Km  are  equal  to  unity,  in 
the  plasma  medium  we  may  write 

Keo  = £1  ' Vo  = yl 


We  shall  now  find  expressions  for  the  azimuthal  component  of  the 
fields  in  terms  of  the  z components  in  the  s'  frame,  i.e.. 


<IV-8>  = %(EZ  ' V 

(iv. 9)  h'  - h'(e;  , 


Assume  the  harmonic  time  dependence  of  the  fields  to  have  the  form 

exp (-ito ' t " ) . Eqs  (IV. 4)  and  (IV. 5)  are  written  in  cylindrical 
coordinates,  and  after  some  substitutions  among  the  component 
equations  of  eqs  (IV. 5)  and  (IV. 5)  we  arrive  at 


(IV. 10) 


U £ / 

O O / 

k'jm  r 


d‘E 


az 
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(VII. 11) 


K'  „ = /i  s2h^  _ a2ay\ 

k_  us  ' 'dz'3$'  3z  ‘ / io  p 


The  results  of  eqs  (IV. 10)  and  (IV. 11)  are  derived  explicitly  in 
appendix  D.  We  may  further  approximate  these  equations  if  we  make 
some  assumptions  about  the  plasma  medium  itself.  First  we  assume 


(IV. 12) 


Es  v Kexp (ikQsin&z) f ($'r 


where  K is  some  constant  that  must  be  determined.  Substituting  eq 
(IV. 12)  into  the  right  hand  side  of  eq  (IV. 10)  yields  a term  that  is 

just  E<j  times  the  factor  yoeQsin  9 If  we  approach  the  cylin- 

der surface  from  the  inside  r-»-a-  , this  factor  may  be  typically  of 
the  order  of  10~1,  for  plasmas  that  are  of  interest  to  us.  However 
an  approach  from  just  outside  the  cylinder  r-a+  produces  a value  of 
sin^Q ' for  the  factor,  and  eq.  (IV. 10)  becomes 


(IV. 13) 
and 

(IV. 14) 


E , = sec 

C 


k ' 2 r 3z'&?' 
o 


nu  cQ  dr 


, r -*-a+ 


E ' « 1 32Et' 

4>  ~ I Ok ' 2r'  JTdf 

o 


1 

' T'1 


r + a_ 


Since  the  azimuthal  field  components  must  be  continuous  across  the 

plasma  surface  we  see  eq  (IV. 13)  and  (IV. 14)  will  be  consistent  if  we 

choose  the  proper  coefficients  A',  B'  >C'  and  D' . 

^ n'  n n n 

Consequently  Eqs  (IV. 10)  and  (IV. 11)  may  then  be  written 

(IV. 15)  E;  = + — — jH  z' 

kQysr'  dz'd<p'  iu't  3r' 

■ iv  16)  H ' = U°e°  + 

$ ^ 2 ^ iu'£  o"r 

using  the  boundary  conditions,  eq  (IV. 1)  and  eqs  (III. 2),  (III. 3), 
and  (III. 10)  and  (III. 11)  substituted  into  eqs  (IV. 15)  and  (IV. 16) 
we  get 
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(IV. 17) 


nsinS'j  (k'acos&0  = - nsinO'ri^  (k"acoso”)A^  + 
no  non 


+ k'  acos6  " H'(1)//  ( k "acos6  ')  B ' + n yS  cos2e  'sin9  'J  ( A 'a)  C ' 
o n o n a n n 


"(t2)  cosV(r)A'aJn(  '' 


A "a)  D' 


(IV. 18)  -k'acosO'J' (k'acose")  = k'acose 'H ^ ^ (k 'acos6 ' ) A' 

o no  o no  n 


- nsin0 'H ^ (k'acosO ' )B 
n o 


; - ($’ 


~ cos29>(\"a) J^(A'a)C^  + 

^o 


+ nsi„e-(^)2< 


cos“6"J  (A'a)D' 
n n 


where  we  have  used  the  approximation, 

(IV.  19;  ~[J)  cos2Q' 

Eqs  (IV. 17)  and  (IV. 18)  together  with  eqs  (IV. 2)  and  (IV. 3)  form  a 
system  of  equations  that  completely  specify  the  undetermined  coef- 
ficients. The  system  of  equations  may  be  written  as  a matrix  equation 

(IV. 20)  AA  = 3 


where 
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We  may  solve  eq  (IV. 20)  for  the  column  vector  A by  left  multiplying 
eq  (IV. 20)  by  a“l,  the  inverse  of  a which  is  defined 


(IV. 21) 


det  | a 


all  a21 ' 


The  factor  det  |a|  is  just  the  determinant  value  of  the  matrix  a.,  the 
primed  elements  a{j  are  just  the  transpose  of  the  cofactors  of  the 
elements  associated  with  a.  Specifically  the  elements  are 

a..  = (k  acosS") 

11  n o 

a12  = 0 


a13  = -Jn(ra) 


a14  = 0 


a21  = 0 

a__  = H ^ (k'acosQ') 
2 2 n o 


-Jn(A'a) 


-nsinO'H^  (k'acosO') 
n o 

k'acosO'Hl^ (k'acosS' 
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a33 
a34  = 
a41  = 
a4  2 = 
a4  3 = 
a44 

and  for  the 


(k '/X " ) 2cos“0 'nJ  (X'a)sinO' 
o n 

- (k'/X')  2cos20'  (X'a)  J '(X'a)  (y./u  ) 

o n 1 o 

k'acosS 'H  ^ ^ (k'acosO') 
o n o 

- sin9'  nH  ^ (k'acosS  ' ) 

n o 

- (k^/X')2cos20'(e1/eo)  ( X'a)  Jn'(  X 'a) 
(k^/X ") 2cos“0 "sin0  'nJn ( X 'a) 

elements  of  the  column  vector  8 , 

- J (k'acosG") 

n o 

0 

sinO'nJ  (k'acosO') 
n o 


|i 

1 


b.  = - k'acosO 'J' (k 'acosO  ') 

4 o no 

The  formal  solution  to  the  scattering  problem  of  a plane  wave  incident 
on  a stationary  plasma  column  is  given  by 

(IV. 22)  A = A-1  8 

It  is  necessary  to  make  a Lorentz  transformation  to  the  observers 
frame  S which  sees  the  plasma  column  translating  with  velocity  ve  . 

One  then  obtains  the  scattered  fields  from  a moving  plasma  column. 


The  resultant  fields 

are 

(IV. 23) 

E(s)  - 

e(!>' 

Z 

z 

(IV. 24) 

H(S)  - 

z 

z 

(IV. 25) 

p(s)  _ 
4> 

V«*S>'-  ''^oKrS-r> 

(IV. 26) 

H<s)  = 
9 

v«#r  + VoEr?,"> 
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( c ) 

IV. 27)  = 

V i.  «<?>') 

r 

r 

z o cp 

(IV. 28)  H(S)  = 

Y(H<!r- 

v e ) 

r 

Z O 0 

These  results  are 

derived  in 

Appendix  E 

V.  THE  DIFFERENTIAL  AND  TOTAL  SCATTERING  CROSS  SECTIONS 

The  differential  scattering  cross  section  is  defined  as 


(V.l) 


da  _ Energy  radiated/unit  time/unit  solid  angle 
dft  " Incident  energy  flux/unit  time/unit  area 


( V.  2 ) 


dSi  = dA/r  = d^dz/r 


is  the  element  of  solid  angle  in  a cylindrical  coordinate  system.  The 
incident  energy  flux  is  the  magnitude  of  the  time-averaged  Poynting 
vector  for  a plane  wave 


(V.3) 


(i)  I _ y 

I — 2? 


4 i Re  lE 


(ilc 


H 


(i) 


)l 


.where  the  symbol  'Pg'  means  tine  real  part  of  the  quantity  that  follows  and  the 
horizontal  bars  denote  the  magnitude  of  the  quantity  between  them. 
For  the  incident  plane  wave  in  our  problem  it  is 


(V.  4) 


s( 


u 


Fr|2Y 


More  specifically,  for  the  case  of  the  infinite  plasma  column 


( V.  5 ) 


d2q 

d<pdz 


limit 
r + « 


w 


|feI2v 


where 


( V.  6 ! 


|§  | = %|Re(E(s)  x H 


(s) 


represents  the  energy  radiated  in  the  radial  direction  out  from  the 
cylinder.  The  component  of  the  Poynting  vector  in  the  azimuthal 
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direction,  S • e^,  will  not  contribute  to  any  real  energy  flow,  but 
instead  represents  some  type  of  "inductive"  energy  that  is  not 
radiative.  The  z component  of  the  Poynting  veccor  has  no  meaning  for 
the  problem  since  the  plasma  is  of  infinite  extent  in  the  z direction. 
The  remainder  of  this  chapter  will  be  spent  deriving  an  explicit  ex- 
pression for  the  scattered  energy  flux  as  given  by  eq  (V.6)  using 
eqs  (IV. 23)  through  (IV.  28)  we  find 

( V.  7)  (E(s)x  H*(s))r  = y(H*ls)  - E*®5 

- w (p  H*<s)  V!r+  e E(f}  V<sr) 
z ^ o z r o z r ' 


(V.  8) 


,-*-(  s)  ->*(s) 
(Ev  'x  H 


- h*!s)  'e<?> 


- YV 


(e  E(s)V<s)%  u H*(s)  V?1') 

Z ' O Z d)  n di  > 


(V.9)  (E(s)x  H*(S))Z  = Y2(E^}  H*<s)  - E<!}  H*1S)  ) 

+ Y2vz(uo[|H^!)  1 2+  |H^r|2]  + eo[|E^}  |2  +|E^}  |2]) 

+ ^zVo  (H^rEMs)  - Hr<s)'EJiS)) 

If  we  use  eqs  (IV. 15),  (IV. 16),  (D.15)  and  (D.16)  and  substitute 

eqs  (III. 10)  and  (III. 11)  into  them  we  find 


(V.10)  (e(s)x  H*(s))r  = 


- i 2/  o 


/ V 

o'  \ o 


X 1 + 


vzsin<0  \^L 


\yr_l,n+n'B^B*;ei(n-n')  (♦'-ir/2)  + 


n.n  =-» 


\ c / n n 

n n ' = - ■» 
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(V.ll) 


(V. 12 ) 
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(E(S)X  H*(s))  = Y 


V 

It,.i2/cc\2  /4isin4>'\ 

; FE  1 ( t ) Uk  r ) 

\ o/  \ o / 


,,n+n\'D*'  i(n-n  ) (<?  -7T/2) 


x AnBn'e 

n ,n'=-30 


l*1"  *»*<•>, 

^sinS'  [i  "(■#)]  + -#  t1  * sin%]) 


n+n'rR-  *'  ^ ^-n')  (<p'-ir/2) 


S«-l'n+n  KV  + BnBn;) 


n,n  =-*> 


where  we  have  used  the  asymptotic  approximation  to  the  Hankel  function 
in  eqs  (V.10)  through  (V.12),  i.e. 


( V.  13)  H JX)  (X) 


exP{i 


{i[;z  - ^ (n+  V)  ] } 


for  large  X,  valid  for  the  far  field  scattering  zone. 


VI.  THE  SHORT  WAVELENGTH  LIMIT  WITH  THE  PLASMA  COLUMN  STATIONARY 
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which  is  of  concern  to  us.  The  results  of  eqs  (VI. 1)  and  (VI. 2)  are 
found  from  eqs  (V.5)  and  (V.10).  For  an  incident  plane  wave  with  a 
frequency 

U>  = 163  GHz 

o 

a plasma  column  radius 
a = . 3 cm 

with  plasma  index  of  refraction 
n = i 3.45 

and  a plasma  frequency 

oj  = 56  3 Ghz 

P 

we  find  that  the  series  in  the  right  hand  side  of  eq  (VI. 2)  converges 
quite  rapidly  after  two  terms,  namely 


(VI. 3)  A'  ~ A*.  * B'  = B* . 

n n n n 


0 , for  | n j , | n' | > 1 


We  also  note  the  relations 


(VI. 4) 


A ' = A' 

-n  n 


B ' = -B' 

-n  n 


We  find  for  normal  incidence  6 = 0C 


(VI. 5)  (2bs>90  = 0U  = 5.03  cm 


A general  expression  may  be  found  for  the  differential 
scattering  cross  section  by  using  eq  (VI. 1)  along  with  the  approxi- 
mation of  eq  (VI. 3)  giving 
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(VI. 7) 


d2a 

2 ( , 2 . | 2 , ,2  f 

cos2  ( $ -tt/2 )] 

d ; ds 

= — < A + B ' “ + 2 i A-  I M 1 + 

mKq  | o 1 1 o ‘ X ■ 

+ 2 1 3 1 1 

~ 1 - COS2U-7T/2)  -2Re  fh  A*+A,A*] 

J i o 1 1 0| 

| cos  ( (f  — 7T/2  ) 

1 

1 

K> 

W 

(D 

,3oB:+B1Bo  l cos  (cp— TT/2)j- 

For  the  grazing  angle  incidence  the  magnitudes  of  the  coefficients 

An  and  dominate  over  all  others  and  they  are  approximarely  equal, 

with  the  angular  dependence  just  that  of  cos  ( 4;-it/2 ) . The  scattering  is 

maximum  in  the  forward  transmitted  direction,  i.e.  at  p = 3ir/2,  and  a minimum  in 

in  the  forward,  reflected  direction  at <}>  = v /2.  For  broadside  incidence  the 

magnitude  of  the  coefficient  A0  is  much  greater  than  all  others  and 

as  before,  the  angular  dependence  is  cos(g-n/2)  and  qualitatively 

the  scattering  angular  distribution  is  identical  to  the  case  of 

grazing  incidence.  Table  1 lists  the  real  and  imaginary  parts  of 

the  complex  coefficients  A„ , A. , B and  B.  for  grazing  and  broadside 
. j o l o 1 

incidence . 

The  terms  "forward"  scattered  and  "back"  scattered  for  non- 
r.ormal  incidence  should  be  physically  interpreted  carefully.  From 
Snell's  law,  and  conservation  of  momentum  it  is  clear  that  the 
forward  scattered  energy  is  just  that  energy  that  is  reflected  from 
the  backside  of  the  cylinder,  a +•  »n  angle  that  is  identically  equal 
to  the  angle  of  incidence.  The  back  scattered  energy  is  that  which 
is  reflected  from  the  front  side  of  the  cylinder,  at  an  angle 
identically  equal  to  the  angle  of  incidence.  Thus  in  both  cases, 
for  non-normal  incidence,  no  energy  would  be  reflected  back  to  the 
source  point  of  the  incident  radiation.  Of  course  for  normal  in- 
cidence the  back  scattered  energy  would  reach  the  initial  source 
point  of  the  incident  radiation. 

VII.  CONCLUSIONS 

The  velocity-dependent  differential  scattering  cross  section  has 
been  found.  For  the  special  case  of  the  plasma  column  at  rest,  and 
the  incident  plane  wave  frequency  sufficiently  high,  we  obtained 
explicit  expressions  for  the  scattered  or  radiated  energy  tlux.  The 
angular  distribution  of  the  scattered  energy  is  transmitted  forward 
for  both  cases  of  incidence,  i.e.  broadside  and  grazing.  This  is 
not  surprising  for  the  situation  of  grazing  incidence,  for  the 
effective  "k0a”  is  large  enough  that  geometric  optics  effects  begin 
to  show,  namely  the  treatment  of  the  waves  as  rays.  In  such  a case 
one  would  expect  the  radiated  energy  to  be  forward  scattered.  This 
is  sometimes  referred  to  as  specular  scattering.  For  broadside 
incidence  the  effective  kQa  is  zero.  One  would  expect,  as  is  the 


14 

case  for  the  scattering  from  the  perfectly  conducting  cylinder 
most  of  scattered  energy  is  in  the  backscattered  direction. 

However,  since  the  dielectric  constant  for  the  plasma  column  is 
0(1)  we  should  not  expect  to  be  able  to  extrapolate  that  well  from 
the  qualitative  behaviour  of  a medium  whose  conductivity  approaches 
infinity.  The  effect  of  the  strength  of  the  dielectric  constant  on 
the  angular  distribution  of  the  scattered  energy  is  isotropic  to 
zeroth  order  in  the  expansion  coefficients  and  approximately  so 
to  the  first  order.  These  results  compare  favorably  with  those  of 
Yeh^O  except  he  presents  the  angular  dependence  of  the  square  of 
the  z component  of  the  E and  H scattered  fields  to  the  incident 
field  instead  of  the  scattering  cross  section  explicitly.  His 
values  for  koa  and  kQa  are  comparable,  however. 


10.  C.  Yeh,  Op.  cit.  p.  5066 

14.  W.  K.  H.  Panofsky  and  M.  Phillips,  "Classical  Electricity  and 
Magnetism,"  (Addison-Wesley  Publishing  Company,  Inc.,  Readinc, 
Mass. , 1962) , p. 330 . 
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APPENDIX  A.  THE  PLANE  WAVE  EXPANSION 

A Four ier-Bessel  expansion  for  a plane  wave  is  developed  from 
the  expansion^ 


(A.  1)  exp{z  (t'--l)/2t}  =£tnJn(Z) 


With  the  substitutions 


t = - exp(iQ') 


z = k' r'cosO 
o 


;e  establish  the  identity 


(A. 2)  ex 


CO 

xp(-ik'cose'y')  = ^ (-l)nexp{in9'}Jn(k'r'cos6') 


iich  is  just  the  result  we  presented  in  eq  (III. 2) 


15  P.  M.  Morse  and  H.  reshuach,  "Methods  of  Theoretical  Physics’ 
Vol.  1,  (McGraw-Hill  Book  Company  Inc.,  New  York,  i - 
dd  619-620 . 
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APPENDIX  B.  THE  PLANE  WAVE  TRANSFORMATIONS 

The  phase  of  a plane  wave  is  a Lorentz  invariant  quantity  since 
one  associates  phase  with  the  number  of  wave  peaks  or  valleys  (or 
any  other  easily  identifiable  reference  points)  that  pass  a given 
reference  point  in  some  given  time.  Imagine  an  observer  in  the  S 
frame  possesses  a phase  meter,  located  at  point  P(z),  as  shown  in 

Fig.  2.  At  time  z = 0,  the  meter  reads  k • z worth  of  phase.  At 
time  T = t later,  the  wave  has  propagated  to  the  left  and  the  meter 

->  -V 

reads  a phase  of  k'Z-tot.  The  number  of  crests  the  meter  has  counted 
is  given  by 

The 

. . Number 

(S.l)  Qf 

Crests 
Counted 

An  observer  in  another  reference  frame  S'  moving  with  velocity 

V = Ve  is  situated  at  point  p'(z').  At  z = 0 , the  two  reference 
frames  % and  S'  coincide.  At  time  t = t'  the  points  P(z)  and 
P'(z')  are  coincident  and  the  number  of  crests  counted  by  the 
observers  in  either  frame  should  be  identical,  i.e. 


phase 


phase  per  crest 


= (k- z-u»t)/2* 


(B.2)  k-z-ojt  = k'*z'-u)'t' 

This  is  shown  in  Fig.  3. 

Recall  the  Lorentz  transformations  between  two  frames  that  are 
in  relative  motion  along  the  z axis, 


x ' = x 


(B . 3) 


y'=  y 

z'=  y(z-vt) 
t'=  y(t-vz/c2) 


Substitution  of  eq  (B.3)  into  eq  (B.2),  collecting  the  coefficients 
of  x,  y,  z and  t and  then  equating  them  to  zero  (since  x,  y,  z,  t 
are  linerly  independent  variables)  yield  the  following  results 
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( B . 4 ) 


k 


x 


k 

y 


O) 


k' 

x 


Y(k'+JVz/c2) 
Y (u)  '+k' ) 


The  inverse  relations  are 


/ k ' as 

v(k  -uv  /c2) 

Z 

z z 

l(i)'  = 

Y(o-kzVz) 

) k ' = 

k 

X 

X 

ii 

k 

y 

If  c and  9 ' are  the  angles  k and  k'  make  with  the  z axis,  and 

for  light  waves  |k  j ■=  -^/c,  |k>|  = w'/c,  we  may  then  recast  the  firs 
two  equations  of  (B.5)  into  the  form 


(co  ' = yw  ( 1-v^cos  9/c) 

(tan  9'  = sin 9/y (cosQ-v^/c) 


where 


(B.  7) 


kz  = |k  | cos  9 
k ' = |k'  i cos  9 

i 2 


The  first  equation  of  (B.6)  represents  the  Doppler  shift  with  the 
Y factor  an  additional  effect.  Even  at  right  angles,  6=^/2,  there 
is  a shift  in  frequency  known  as  the  transverse  Doppler  shift. 
Relating  these  results  to  those  of  (III. 4)  through  (III. 9)  we  must 
first  identify 


eo  = V2  -e 

The  first  equation  in  (B.6)  becomes 
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(B.  8) 
Now 


u ' = Yu  (1-  gsin9Q) 


k x + k y = k'x'+k'y' 
x y x y 


This  implies 

ksin9  = k'sin6' 

and  for  the  case  of  the  cylinder 
( B . 9 ) k^cosQ'  = kQCos0o 

The  second  equation  in  (B.6)  becomes 

(— JVcosS'  = (cosG-vz/c)  1 

With  some  manipulation  and  the  use  of  the  first  equation  in  (B.6) 
arrive  at 


(B.10)  sine'  = (sineQ-e ) / (l-6sin6o) 

We  have  already  verified  eq  (III. 7)  since  it  is  no  more  than  the 
Doppler  shift  equation  of  (B.6). 

We  must  now  examine  the  field  transformation  equations13 


(B.ll) 


E = E 
II  II 

E'  = Y(E  + v x B) 
1.  X 


More  specifically 


(B. 12) 


E cos  6 e 
o o z 


(B . 13 ) 


= Y(E 


v x B)  = Y(-EQsin0oe  + v x B) 


13.  J.  D.  Jackson,  Op.  Cit.  P223-226. 


we 
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and  consequently 


(B.15) 


v x B = — E e 

c o y 


Substitution  of  eq  (B.15)  into  eq  (B.13)  yields 


(B.16)  E = YEo(3  - sin9Q)ey 

and 

(B . 17 ) E'  2 = E2 (cos20  + y 2sin2S  + Y262)-  2E2y2sin9 

O O o o o o 

With  some  rearrangement  and  explicit  use  of  the  defining  equation 
for  y we  obtain 

(B . 18)  E^  = Eqy(1  - 3sin0o) 

Thus  eq  (III. 8)  is  confirmed. 

The  equality  of  phases,  eq  (B.2)  combined  with  the  first  equation 
in  (B.ll)  gives  rise  to  eq  (III. 9).  We  have  now  established  eqs 
(III. 4)  through  (III. 9)  of  Chapter  III. 


26 


NSWC/WOL/TR  76-101 


APPENDIX  C.  THE  WAVE  NUMBER  IN  THE  PLASMA 

If  we  consider  the  incident  plane  wave  striking  a tangent  plane 
of  the  plasma's  cylindrical  surface  we  immediately  obtain  Snell's 
lawi? 


(C.l) 


>in0 VsinG " = / Wo2,. 


where  0"  and  0"  are  the  angles  the  propagation  vectors  and  kg' 

make  with  the  normal  to  the  tangent  plane  outside  and  inside  the 
cylinder  respectively.  Define 


(C.2)  kg  = (k^'  cos0")/cos0"  = A'/cos0" 

Since  the  S'  frame  is  identically  the  same  inside  or  outside  the 
cylinder 

(C . 3)  w't'  = u^t" 

where 

(C . 4 ) w'  = k'c 

o 


(C.5)  til"  = kg'c  = Ad/cos 0" 

Since  t'  is  identically  equal  to  t"  we  obtain 

(C.  6) 


■%  _ 


“T  = k^/(w0e0)  2 = A '/(p1£1)2cos6" 


or 


a-)'  ■(£) 


COS^0" 


We  may  rearrange  eq  (C.l)  into  a more  convenient  form  by  squaring 
it  and  using  a tregonometric  identity  to  obtain 


(C . 8) 


V-oW 


2a. 

!COS  0 


ylel 

yo£o 


- sin^  0 ' 


Substituting  the  left  hand  side  of  eq  (C.7)  into  eq  (C.8)  gives 
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APPENDIX  D.  THE  TRANSVERSE  FIELD  COMPONENTS  AS 
FUNCTIONS  OF  THE  LONGITUDINAL  FIELDS 

Assuming  harmonic  time  dependence  of  the  form  exp(-i  't')  the 
Maxwell  equations  become 

(D.l)  7'x  E'=  &'B'  = iu"<mu  u'  , V'-H'  = 0 

(D. 2)  7"x  H'  = -i-j'<ecE  ' , V • E'=  0 

For  any  vector  field  A (r)  we  may  write  its  curl  in  cylindrical 


coordinates 

as 

-> 

1 3AZ 

3 A 

(D.  3) 

(V 

X 

A)  r = 

r 3<J> 

3 r 

3Ar 

3AZ 

( D . 4 ) 

(7 

X 

A)  = 
<P 

3 z 

3r 

(D . 5) 

X 

II 

N 

1 3 
r SF 

1 3Ar 
r 

then  the  first  equations  in  (D.l)  and 


(D . 6) 

i'j  Kmu0H^  ; 

1 3EZ ' 
r ' d<p  ' 

3E* 

~Jz 

(D.  7) 

dE?'  3E2 

i u'<mu„H  ' , = 
o <p 

3z  ' " TT 

(D.  8) 

k ^ (r' 
r 3r 

E%> 

0 

and 

( D.  9 ) 

1 3HZ/  3 

HI' 

c'  do' 

3 z ' 

(D.10) 

-iu)'<£QE'  + 

3H2  - _ 3HZ, 
3z'  ” 3 r ' 

(D . 2 ) 


3Er 

oO 


become 
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1 a i U11r 

(d. id  - i-  ■£=-  «*'»;-)  - — 


Since  we  are  first  attempting  to  establish  the  relations 

Ea  = Ea  (E  , H ) , H , = H . (E  , H ) 

4>  z z <j>  <J)  z z 

We  shall  differentiate  eq  (D.6),  then  substitute  this  result  into 
eq  (D.10) . We  find 


ye  ~ 

(D.  12)  E'  - -2-2  5 

t i.  ^2.  - v r c 


- 2 E - 
o z 


V ^ ' 


Tz^SV 


V'\  , 1 3H\ 

V 1U>  e — 


where  we  have  used  the  relation 


(D.13)  u>'2  = k'2/yoeo 


Let  us  differentiate  eq  (D.9)  and  substitute  the 
result  into  eq  (D.7)  getting 


(D. 14 ) H 


3eo  (l  _ 92h:A  _l 

^"ue  ^ r ' ?V3b'  T7r'^—j  ' i - 


. y Sr 


In  a similar  fashion  one  may  find  the  radial  components  of  the  fields 
in  terms  of  the  longitudinal  components,  i.e. 


Er  = VW  ' Hr  Hr  ''Ez,Ez) 

Differentiate  eq  (D.7)  with  respect  to  z'  and  substitute  the  result 
into  eq  (D.9).  We  obtain 


(D. 15)  E , 


r m e l r 


3HZ-  po£o  33Hz'  

34>'  KqVe r ' 3‘z'3<T 


32EZ'  \ 
"3  z ' 3 r j 


and  we  have  established  the  first  relation. 


Differentiate  eq  (D.10)  and  substitute  this  result  into  eq  (D.6). 
The  result  is 


30 


o o 


NSWC/WOL/TR  76-101 


APPENDIX  E.  THE  SCATTERED  FIELD  COMPONENTS 
IN  THE  OBSERVERS  FRAME 

If  we  use  the  Lorentz_ transformation  equations  appropriate  to 
the  electromagnetic  field-^--’  we  see 


(E.l) 


E = E ' 

I!  II 


B = B ' 
II  II 


E = E'  . 
z z 


B = B'  . 
z z 


(E . 2 ) Ej_  = y(E±  - v x B±  ) 


For  the  case 


E.  = E , e.* 

1 4>  <P 


<v  * MV  • W5< 


and  we  have 


<E-3>  E*  = VE'.  - V2U0Hr) 


For  the  other  transverse  component 


E,  , = EVe  . 
i r r 


( x Bi.)r.=  - vzB'.er, 


which  gives 

(E • 4 ) Er  = y(E^  + vzuoHj.) 


W.  K.  H.  Panofsky  arid  M.  Phillips,  "Classical  Electricity  and  Magnetism," 
(Addison-Wesley  Publishing  Company,  Inc.  Reading,  Mass. , 1962) , p.  330. 
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